Evidence of phonon mediated coupling between charge transfer 
and ligand field excitons in Sr2Cu02Cl2 



R. Lovenich, A. B. Schumacher!, J. S. Dodgei, and D. S. Chemla 

Department of Physics, University of California at Berkeley, Berkeley, California 94720 
Materials Sciences Division, E. 0. Lawrence Berkeley National Laboratory, Berkeley, California 

94720 

L. L. Miller 

Ames Laboratory and Department of Physics, Iowa State University, Ames, IA 50011 

(February 1, 2008) 



Abstract 

We present a comparative experimental and theoretical investigation of 
the two-dimensional charge-transfer gap in the strongly correlated material 
Sr2Cu02Cl2. We observe an Urbach behaviour in the absorption profile over 
a surprisingly wide range of energies and temperatures. We present a model 
that accounts for phonon scattering to infinite order and which allows us to 
explain the data accurately by assuming coupling of the charge transfer gap 
exciton to lower energy electronic excitations. 



I. INTRODUCTION 



The excitation spectrum of transition metal oxides is characterized by strong coupling 
of spin, charge, orbital and lattice degrees of freedom, and the complexity of these interac- 
tions has severely hindered progress in understanding the phenomenon of high-temperature 
superconductivity in the layered copper oxides. Even the basic energy level spectrum of the 
undoped, insulating cuprates continues to be controversial. The dominant optical absorption 
feature in these materials is a broad excitonic peak at the charge-transfer (CT) gap, about 
half an electron volt wide, which exhibits strong dependence on temperature and dopingii. 
Despite the strong interest in insulating cuprates and the fundamental physical significance 
of the CT gap exciton, few attempts have been made to explain its structure and coupling 
to phononsii. 

In this paper, we present an experimental investigation of the linear absorption of 
Sr 2 Cu02Cl2 and a new theory for the CT excitation in the CuO-plane. The theory in- 
cludes the coupling of the exciton to phonons and to an additional, low energy electronic 
continuum of states. It explains the linear absorption of Sr2Cu02Cl2, and in particular the 
Urbach behaviour over a wide range of temperatures and energies. 



II. SAMPLE STRUCTURE 

Sr 2 Cu02Cl 2 is a two-dimensional, spin-1/2 Heisenberg antiferromagnet with a Neel tem- 
perature Tjy ~ 250 K and an exchange interaction energy J=125 meV0. In the ground state 
a single hole per Cu02 unit cell is located on the Cu site and exhibits a 3d x 2_ y 2 symmetry. 

The lowest energy optically allowed electronic excitation in this material corresponds to 
the Cu 3d x 2_ y 2 to O 2pa transfer of a hole in the Cu0 2 layer. This transition classifies the 
material as a charge-transfer insulator!!. The optical absorption spectrum consists mainly 
of a rather sharp peak near the band edge at about 2 eV and another wider peak at about 
2.5 evi A detailed analysis of the energy region below the band edge reveals a weak ab- 
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sorption step at about 1.4 eV, two orders of magnitude smaller than the main CT structure. 
It is attributed to a Cu ligand field transition!, which should be optically forbidden, but 
acquires a small oscillator strength owing to a breaking of the crystal symmetry. This transi- 
tion has also been observed in resonant X-ray Raman spectroscopyi and large-shift Raman 
scattering!. 

The dominant peak at the band edge is attributed to a bound excitonic state. A model 
by Zhang and Ngi predicts an exciton formed by an electron at a Cu site and a delocalized 
hole at the four neighbouring O sites, with a stronger weight at the O site from which the 
electron came. The model gives a large bandwidth for the center of mass motion of this 
exciton of about 1.5 eV0. This means that the exciton has an effective mass smaller than the 
electron or hole alone, which can be understood by the fact that the transfer of the exciton 
from one site to another does not disturb the antiferromagnetic order of the background. 



III. EXPERIMENT 

We used high quality single crystals of Sr2Cu02Cl2 which were grown by the technique 
described in Ref. |ll]. This method yields platelets with a fairly large (001) surface. By 
tape-cleaving the crystals, we reduced the sample thickness, d, to below 100 nm. These thin 
samples enabled us to measure the absorption coefficient, a(u), across the charge transfer gap 
directly, by comparing as a function the frequency u the intensity transmitted through the 
sample on the substrate, Itm, with the intensity transmitted through the substrate only, Isub- 
a(uj) = ~2 ^ n (jf^)- This procedure is more precise and direct than the previous reflectivity 
measurements, which relied on Kramers-Kronig transformation to extract information about 
the absorption coefficient!^!. Figure [I] shows the linear absorption coefficient for selected 
temperatures. The crystal used for this measurement was d = 95±5 nm thick and held 
inside a cold finger cryostat. The inset details the weak onset of the absorption at 1.4 eV 
and was measured using a 300 /im thick sample. 

The excitonic peak shows a strong temperature dependence. With increasing tempera- 
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ture the width — measured as the half width at half maximum — increases from 180 meV 
at 15 K to 270 meV at 350 K. The maximum exhibits a redshift of about 80 meV in this 
temperature range. The temperature dependence of the excitonic peak is shown in figure |2|. 
Both the shift towards lower energies and the broadening can be fit by a Bose-Einstein 
occupation function with a single oscillator!. From this fit we find the energy of this oscilla- 
tor, huo ~ 45 meV, to agree with that of a previously identified bond-bending longitudinal 
optical (LO) phonon0 and well below the characteristic magnetic energy of J=125 meV 
in Sr2Cu02Cl2. Experimental evidence therefore suggests coupling of the CT-exciton to 
phonons rather than to magnetic excitations as the origin of its red shift and broadening!. 
This is further supported by a careful analysis of the line shape of the low energy side of the 
CT excitation. As shown in figure ^ a(uS) exhibits an exponential behaviour over a wide 
range of temperatures and energies. We find the experimental data can be fit very well to 
the Urbach formula 

a(T)(E-E ) 

a(u,T) = a exp — ; (1) 

K B 1 

where Eo and ao are temperature independent parameters, and Tiujq is the energy of the 
phonon to which the exciton couples. The parameter <7o yields information on the strength 
of this coupling and Eo is related to the strength of the exciton binding. From a global best 
fit to the data plotted in figure |3] data we obtain fru;o=45 meV, E =1.95 eV and <7o = 0.35. 
The latter indicates a very strong coupling between excitons and phonons. Furthermore, 
the fact that Eo is smaller than the energy E(a=max) where the absorption is maximal, 
suggests a tightly bound exciton0. 

It is worth noting that it is very surprising to observe an Urbach behavior at temperature 
as low as 15 K where the occupation probability of LO-phonons with ^a; =45 meV is 
Nq f«2xl0 -16 and absorption of thermal phonons is negligible. 
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IV. THEORY 

To the best of our knowledge, the model of Falck et al.l is the only one in the literature 
describing the temperature dependence of the CT excitations. Falck et al.l investigated the 
absorption line shape of a related compound, the charge-transfer insulator L^CuCU. The 
excitonic peak in this material also shows a temperature shift and broadening. Approximat- 
ing the Coulomb interaction with an attractive contact interaction between the electron and 
hole, they derived a simple expression for the absorption coefficient. While without interac- 
tions the interband contribution to the absorption is just the density of states, the Coulomb 
interaction leads to an increase of the absorption near the band edge. This describes the 
basic shape of the spectrum. Falck et al. account for the shift of the band edge and the 
finite lifetime of the exciton by inserting the real and imaginary part of the one-phonon self 
energy. With four parameters they are able to fit well the linear absorption in L^CuO^ 

Our experimental results show the exponential behavior of an Urbach tail down to very 
low temperatures. It is in striking contradiction with the prediction of Ref. 0, which by 
construction always produces a Lorentzian tail. The coupling constant used in Ref. |2| was 
a = 10. 8^. This is an extremely high value compared to normal semiconductors (where it is 
well below 1.0) and even other materials with a strong electron-phonon coupling. For such 
a high coupling constant the one-phonon approximation is no longer valid. 

Since a strong coupling to LO-phonons seems to be the reason for the strong temper- 
ature dependence of the spectra, we propose a more advanced model, taking into account 
many-phonon processes. The approach was inspired by results obtained for conventional 
semiconductors!!!, where Urbach tails are found for high temperatures. 

The theoretical analysis is organized as follows: First we introduce the Hamiltonian in 
an exciton basis. The Green's function of that Hamiltonian will be calculated by means 
of the cumulant expansion. We will then discuss the dependence of the results on model 
parameters for different cases of physical interest. We will show that it is necessary to 
include ligand field (LF) excitations, which were not yet known at the time of Ref.i. These 
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are the states into which the CT exciton can scatter even at low temperatures by phonon 
emission, thus leading to the Urbach tail. 
The model Hamiltonian: 



H = H + V 



(3) 



#0 = E E nM fl n!k <k + fad E b \ 6 q 



(4) 




(5) 



is written in an exciton basis, where a^ k is the creation operator for an exciton of type 
j, internal quantum number n and center of mass momentum k. In the case of the j=CT 
exciton we account for a bound state (n — 0), and for the unbound exciton electron-hole 
continuum (n > 0). The pair (n, k) can then be transformed to the pair of electron and 
hole momentum, which is more suited to describe the free electron-hole pair. In the case 
of the ligand field excitation (j—LF) we consider continuum states only. The operators 6 q 
and respectively describe the annihilation and creation of an dispersionless LO-phonon 
of frequency u . 

The imaginary part of the dielectric function and hence the spectral features of the 
absorption coefficient are given by the imaginary part of the retarded Green's function 



where the summation includes only the optically active CT-exciton. 

Since the excitation starts from a ground state characterized by empty hole and electron 
states, the retarded Green's function equals the time-ordered one, for which an expansion 
in the electron-phonon coupling exists: 




j=CT>k 



E 




(6) 



• oo 

G>(nk;f) =-^9(t) £W2(nk;f) 



(7) 




' o o 



Jdt 2v (Tai ik (t)V(t 1 )...V(t 2 „)a% k (0)) , 



(8) 
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with the time-ordering operator T. The unperturbed Green's function is given by 
G 3 (nk;t) = —i/h9(t) exp(—iE 3 nk t/h). The first order term in (|7|) would give rise to one- 
phonon processes. However, it is possible to account for multi-phonon processes by resum- 
ming the expansion Eq. (|7]) in form of an exponential, leading to the cumulant expansion or 
linked-cluster theory. We write 

G*(nk; t) = Gi(nk; t) exp ( £ f*(nk; t) ) , (9) 

where F u denotes a contribution which contains the 2z/-th power of the coupling matrix- 
element. Comparing the different powers of the expansions (0) and (|Sp the following relations 
hold: 

F((nk;t) =e iE U t/h Wi(nk;t) (10) 
Fi(nk; t) = e iE U t/h Wi(nk; t) - i F( 2 (nk; t) . (11) 

The advantage of this re-summation is that F\ already includes independent phonon scatter- 
ing to infinite order. Thus as long as they are independent from each other, the first term in 
the cumulant expansion F\ accounts for all these contributions in contrast to the expansion 
(0). Similarly F 2 describes all electron-phonon interactions involving simultaneously two 
phonons, etc. The cumulant expansion (||) is expected to converge much faster than (0) 
and to be appropriate for low and intermediate coupling constants. Up to now we restricted 
ourself to the lowest order only, which means that we treat the fields induced by the lattice 
distortions classically. 

The time-integrations in (R) can be performed analytically, yielding: 



F((nk = 0;t) = £ |M£ iq | 2 {(N + l)/« q - < + hoj ,t) 

h'mq 

+N f(E£ M -Ei i0 -huo,t)}. (12) 



The time dependence is described by the function 



f(s,t) =l(l- e -^-^). (13) 



This function is regular for all values of e and has its maximum value at e = for finite 
times, which in Eq. (0) corresponds to E^ = E 3 n Q ±hcu , i.e. the absorption and emission 
of LO-phonons in the (j', mq) <-> (j, n) transition. At low temperatures only the emission 
term proportional (Nq + 1) — > 1 contributes, which peaks out one phonon energy below the 
CT exciton. 

In order to keep our model simple, we neglect the q dependence of the coupling matrix 
elements as well as the dependence on the relative momentum if m > or n > 0, 

i.e., for the electron-hole continuum. This means we have to consider the following five 
scattering processes: 

M xx Scattering from the CT-exciton into its center of mass continuum; 

M xc Scattering from the CT-exciton into its own electron-hole continuum; 

M cc Scattering within the CT electron-hole continuum; 

M x d Scattering from the CT-exciton into the LF exciton band; 

M c d Scattering from the CT electron-hole continuum into the LF exciton band; 

Since the matrix-elements do not depend on momentum variables any more, the sum- 
mations over momentum variables in flT2|) can be transformed into energy integrals with the 
density of states into which the CT-exciton is scattered: 

F 1 (n = {x,c},0;t) = Jde ]T M nm z m {e) (14) 

m={x,c,d} 

x {(N + l)/(e - E nfi + huo, t) + N f{e - E nfi - hu , t) } . 
The densities of states are defined as: 

z ^ £ ) = -^Y, 6 ( £ ~ E x,a) and ( 15 ) 

q 

4E^- £ m). (16) 

qk 
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We model the energy dispersion of the exciton by a two dimensional tight binding dis- 
persion E = B/4 - (2 — cos q x — cos q y ) with a bandwidth of 5=1.2 eV and a electron-hole 
relative motion continuum, also with a 2d tight binding dispersion, but a 0.6 eV bandwidth. 
We have solved numerically Eq. (|6|) with this dispersion law in three cases of physical in- 
terest. In case (1) we consider the situation where phonon scattering occurs only within the 
CT-exciton states, i.e., CT bound-exciton and its electron-hole continuum, and the transi- 
tions into other exciton bands are not allowed, M x d = and M c d = 0. In the other cases 
we assume that phonon scattering can also occur with other excitonic states0! (LF exciton 
band) of energy lower than the CT exciton, M x d ^ and M c d ^ 0, and to explore the 
generality of our assumptions we consider two slightly different LF exciton bandwidths in 
the cases (2a) and (2b). 

The results in case (1), for the contribution of the bound CT-exciton (n = x in @) to 
the absorption only, are shown in figure f|and are in obvious contradiction with experiment. 
Clearly the coupling to LO-phonons leads to a broadening of the CT-exciton line. It be- 
comes asymmetric with a longer tail at the high energy side. On the low energy side, the 
spectrum has an Urbach tail, but the width is extremely small, especially at low tempera- 
tures. This is true even for very high coupling constants. Accounting for the large width at 
low temperatures by assuming an homogenous or inhomogenous broadening would lead to 
a Lorentzian or Gaussian profile of the lower energy side of the spectrum. 

The shortcoming of our model in case (1) is clear. Since the Urbach behavior is observed 
at very low temperature phonon emission governs the interaction with the lattice as implied 



by Eq. (0). Therefore, for broad CT-exciton peak with a flat exponential tail to occur 
there must be states available below its energy and the CT-exciton must be able to couple 
to these states by phonons. Only then can LO-phonon scattering be strong even at zero 
temperature. Both our own data and the experiments in Ref. |^ show that such a continuum 
exists and starts at about 1.4 eV. It is assigned to Cu d-d transitions and can couple to the 
excitons by means of LO-phonons. 

Therefore we expanded the model by including the continuum of LF excitations. Since 
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little is known about the bandwidth of these excitationsEDtl we have considered two values 
for this parameter. In case (2a) we assumed that the LF bandwidth is the same as that of 
the CT electron-hole continuum. The CT exciton binding energy was fitted to be 240 meV, 
and \M nm | 2 = 20(huo) 2 was used for all the coupling constants. In addition to the LO 
phonon broadening we included an inhomogeneous broadening T/=60 meV to account for 
sample imperfections. An homogenous broadening also should be included, however, it 
must be very small compared to Fi because otherwise a Lorentzian profile would show up 
for energies far below the CT exciton energy in contradiction with experiment. Therefore, in 
order to limit the number of parameters we have not included it in our calculations. Finally, 
to fit the spectral contribution of the CT continuum we found that the spectral weight of 
the electron-hole continuum needs to be twice that of the exciton. 

The solid lines in figures |5] and |6| show the theoretical results in case (2a) both on a linear 
and logarithmic scale. They are in excellent agreement with the experimental data, shown 
by the dots. The broadening and the shift of the CT exciton are reproduced correctly and 
even the decrease in oscillator strength is well described. 

It may be argued that the LF excitations bandwidth should be smaller than that of the 
CT electron-hole continuum. In order to investigate this eventuality, we have considered 
in case (2b) the effects of a narrower LF bandwidth, only 1 eV, while maintaining the flat 
density of states starting at 1.4 eV in agreement with the experimental absorption spectrum. 
In that case we find that we can also achieve an excellent agreement with experiment, of 
the same quality as in case (2a) (not shown), but then we need two independent coupling 
constants, one describing the phonon coupling to the LF excitations and another coupling 
within the CT manifold. Interestingly, the numerical integration points out that the tem- 
perature shift of the CT exciton is governed by the coupling to its own continuum, whereas 
its broadening is mostly determined by the coupling to the LF excitations. With these two 
parameters the model is also robust against changes of the onset of the LF excitations. 

Although the model formally has a number of fit parameters, many (phonon energy, CT 
exciton bandwidths, onset of LF excitations etc.) are fixed by experiment. Furthermore, we 
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have limited the number of parameters by limiting the number of phonon coupling constants 
to one (case-2a) or two (case-2b), so that in fact only few parameters are free: the (one or 
two) phonon coupling constants, the inhomogenous broadening, the relative spectral weight 
of the CT exciton and its continuum, the CT exciton binding energy and an overall scaling 
factor. 

The results indicate that the coupling constants thus obtained are indeed large. Due 
to the different models they cannot be directly compared to the results of Falck et al. It 
is straight forward to check if accounting only for independent one-phonon processes is 
sufficient, since one can calculate F2 with not too much effort if the q-dependence of the 
matrix elements is neglected. This has not yet been done. Another theory often used to 
describe localized electronic excitations coupled to LO-phonons is Toyozawa's@. In the 
strong coupling regime it predicts a self-trapping of the excitons, i. e. due to the phonon 
cloud surrounding the exciton its effective mass becomes extremely large. Since we have 
shown on very general physical grounds that the low temperature Urbach behavior requires 
the coupling of the CT-exciton to a continuum of states below the exciton energy, we expect 
that a model based on Toyozawa's theory would give results qualitatively similar to ours. 

In conclusion we have presented measurements of the linear absorption of the material 
Sr2Cu02Cl2 which exhibits broadening and shift of the charge transfer gap exciton and, 
importantly, an Urbach behavior down to very low temperatures. This behavior cannot 
be explained by the current theories. We have developed a model based on a cumulant 
expansion that explains well all our data providing that the charge transfer gap exciton 
is coupled by LO-phonons to a continuum of states with lower energy. We identify this 
continuum with the Cu d-d continuum which is revealed in the linear absorption spectra 
owing to a weak crystal field symmetry breaking. A consequence of our interpretation is that 
the magnetic degrees of freedom do not seem to play a significant role in the temperature 
dependence of the charge transfer absorption edge in this material. 
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FIGURES 
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FIG. 1. Absorption coefficient of Sr2Cu02Cl2 at different temperature, measured on a 
1=950 A thick sample. Inset (l=300^m): the weakly allowed d-d transition. Solid lines T=15 K, 
dashed line T=150 K, dotted lines T=250 K, dashed-dotted line T=350 K 
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FIG. 2. Half width (solid circles, left axis) and energy (open circles, right axis) of the absorption 
peak at different temperatures. The dotted shows the best fit to a single oscillator Bose-Einstein 
occupation function with Efe oson =45 meV 
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FIG. 3. Logarithmic plot of the absorption coefficient of S^CuO^Ckj at selected tempera- 
tures, displaying the Urbach-tail. Solid squares T=15 K, open squares T= 150 K, open triangles 
T=250 K, solid circles T=300 K, open circles T=350 K, diamonds T=400 K. The solid lines are 
fits to the Urbach formula. 
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FIG. 4. Absorption coefficient for the excitonic contribution only on a logarithmic scale, cal- 
culated for scattering into the exciton and e-h continuum only. Solid line T=15 K, dashed line 
T=300 K. The parameters are the same as for the final calculations, except that the coupling 
constants are twice as big and the inhomogenous broadening was only 2 meV. 
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FIG. 5. Absorption coefficient data and theory. Only three temperatures are shown for clarity. 
Squares T=15 K, closed circles T=250 K, open circles T=350 K, solid lines are the theoretical 
results. 
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FIG. 6. Absorption coefficient data and theory, displayed on a semi-logarithmic scale. Only 
three temperatures are shown for clarity. Squares T=15 K, closed circles T=250 K, open circles 
T=350 K, solid lines are the theoretical results. 



20 



